ESTIMATION OF THE ROAD CONDITION UNDER A VEHICLE 
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FIELD OF THE INVENTION 

The present invention relates generally to the estimation of 
5 the road condition under a vehicle and, for example, to 
systems, methods, and computer program products for estimating 
the road condition under a vehicle. 

BACKGROUND OF THE INVENTION 

Modern cars comprise electronic control systems as anti-lock- 
10 braking systems (ABS) , dynamic stability systems, ant i- spin 
systems and traction control systems. Besides these active 
control systems there also exist driver safety information 
systems as road friction indicators and sensor- free tyre 
pressure monitoring systems which present information about the 
15 driving condition to the driver. 

All the above-mentioned systems benefit from the knowledge 
about the road aurface condition under the vehicle. Several 
different techniques are used in the prior art to determine the 
road surface condition under a driving vehicle. One such 
20 technique is based on vertical accelerometers in a suspension 
system of a car. Another technique is based on level meters in 
the fuel tank of the car. Other techniques use special air mass 
flow sensors in the engine control unit. 

The present invention relates to techniques for estimating the 
25 road condition which make use of the signals obtained from 
wheel speed sensors, e.g. the wheel speed sensors of standard 
anti-block braking systems. Using the signals from wheel speed 
sensors of ABS systems (and/or from the vehicle's internal CAN- 
bus) provides an economical way to road surface condition 
30 measurements since these ABS systems belong to the standard 
equipment of the majority of the cars and trucks sold today. 

Such a system which is based on the signals of wheel speed 
sensors is for example disclosed in US-patent 5,566,090 which 
is directed to a method for detecting stretches of bad road 
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directly from the raw data provided by an ABS sensor. The 
method uses the fact that stretches of bad road result in 
strong fluctuations of the wheel speeds of the car. Strong 
wheel speed fluctuations in turn result in large differences 
5 between successive segment times, where the segment time is the 
time the wheel needs to pass through associated angle segments. 
The disclosed method determines a stretch of bad road if the 
difference between successive segment times is greater than a 
pre-set limit value. This simple decision algorithm operates 
10 directly on the raw signals of the wheel speed sensor. The US 
4,837,727 discloses a method which is based on a similar 
decision algorithm. 

EP 0 795 448 A2 discloses a road surface condition detection 
system which comprises a wheel speed sensor for detecting a 

15 wheel speed of at least one wheel to generate a wheel speed 
signal and a control unit which integrates the wheel speed 
signal for a predetermined period of time. The control unit 
determines a rough road surface condition when the integrated 
signal is above a predetermined threshold value and, otherwise, 

20 a normal road surface condition. Before the integration, the 
wheel speed signal is band-pass filtered in the frequency range 
of 10-15 Hz. 

SUMMARY OF THE INVENTION 

A first aspect of the invention is directed to a system for 
25 estimating the ground condition under a driving vehicle. The 
system comprises a wheel speed sensor for sensing a wheel speed 
signal which is indicative of the wheel speed of a vehicle's 
wheel driving over the ground and a first analyser unit coupled 
to said wheel speed sensor. The first analyser unit comprises a 
30 sensor imperfection estimation section which is designed to 
estimate a sensor imperfection signal from the wheel speed 
signal which is indicative of the sensor imperfection of the 
wheel speed sensor, a signal correction section which is 
designed to determine an imperfection-corrected sensor signal 
35 from the wheel speed signal and the sensor imperfection signal, 
and a ground condition estimation section which is designed to 
estimate a first estimation value indicative of the ground 
condition from the imperfection-corrected sensor signal. 
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Another aspect of the invention is directed to a method for 
estimating the ground condition under a driving vehicle, 
comprising the steps of: 
5 - sensing a wheel speed signal by means of a wheel speed 

sensor which is indicative of the wheel speed of a vehicle's 

wheel driving over the ground; and 

estimating a sensor imperfection signal from the wheel speed 
signal which is indicative of the sensor imperfection of the 
10 wheel speed sensor; 

determining an imperfection-corrected sensor signal from the 
wheel speed signal and the sensor imperfection signal; and 
estimating a first estimation value indicative of the ground 
condition from the imperfection-corrected sensor signal. 

15 

A further aspect of the invention is directed to a computer 
program including program code for carrying out a method, when 
executed on a processing system, of estimating the ground 
condition under a driving vehicle, the method comprising the 
20 steps of: 

sensing a wheel speed signal by means of a wheel speed 
sensor which is indicative of the wheel speed of a vehicle's 
wheel driving over the ground; and 

estimating a sensor imperfection signal from the wheel speed 
25 signal which is indicative of the sensor imperfection of the 

wheel speed sensor; 

determining an imperfection-corrected sensor signal from the 
wheel speed signal and the sensor imperfection signal; and 
estimating a first estimation value indicative of the ground 
30 condition from the imperfection-corrected sensor signal. 

Other features are inherent in the methods and systems 

disclosed or will become apparent to those skilled in the art 

from the following detailed description of embodiments and its 
35 accompanying drawings. 
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DESCRIPTION OF THE DRAWINGS 

Embodiments of the invention will now be described, by way of 
example, and with reference to the accompanying drawings, in 
5 which: 

Fig. 1 shows a car having four wheels and driving on a road 

which changes in driving direction from a normal 
surface condition to a rough road condition; 

Fig. 2 schematically shows a wheel speed sensor comprised 

io of a segmented rotary element and a sensor element; 

Fig. 3 shows an exemplary diagram of four wheel speed 

signals obtained from the four wheels of a driving 
vehicle as a function of time; 

Fig. 4 shows a diagram representing a wheel speed signal as 

15 a function of the sample number; 

Fig. 5 shows a block diagram of an embodiment of the system 

for estimating the road condition under the vehicle, 
the embodiment comprising a wheel speed sensor and 
an analyser unit ; 

20 Fig. 6 shows a block diagram of an embodiment of the ground 

condition estimation section which is part of the 
system of Fig. 5; 

Fig. 7 shows a block diagram of a further embodiment of the 

ground condition estimation section which is part of 
25 the system of Fig. 5; 

Fig. 8 shows a block diagram of an embodiment of the 

variance estimation section which is part of the 
ground condition estimation sections of Fig. 6 and 
Fig. 7; 

30 Fig. 9 shows a block diagram of an embodiment of the system 

for estimating the road condition under a vehicle 
which is based on the signals of four different 
wheel speed sensors; 
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Fig. 10 shows a block diagram of an embodiment of a decision 
unit of the system for estimating the road 
condition; 

Fig. 11 shows a block diagram of an embodiment of the system 
5 with two different types of analyser units; 

Fig. 12 shows a block diagram of an embodiment of the second 
type of analyser unit comprising a filter section; 

Fig. 13 shows a block diagram of an embodiment of the system 
for determining the road condition, wherein the 
io system comprises two types of analyser units which 

evaluate the signals from four different wheel speed 
sensors ; 

Fig. 14 shows a block diagram of an alternative embodiment 
to the one of Fig. 13; 

15 Fig. 15 shows several diagrams representing the operation of 

the embodiment of Fig. 13. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 shows a car 1 having four wheels and driving on a road 
which changes its surface condition from a normal surface 
20 condition 2 to a rough surface condition 3. A normal surface 
condition may for example be present when the car 1 is driving 
on an asphaltic road. A rough road condition may for example 
occur on gravel, rough asphalt, rough ice and some types of 
snowy roads. The arrow labelled with v in Fig. 1 indicates the 

25 driving direction of the car 1. The arrow labelled with co 
indicates the wheel rotation which is caused by the forward 
movement of the car 1. In the description of the embodiments 
which follows now, the principles of the invention are 
explained with reference to a car having four wheels. However, 

30 the proposed systems and methods may as well be applied to 
other types of vehicles, as for example trucks, buses and motor 
cycles, having a different number of wheels. 
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Wheel speed sensor imperfections 

Fig. 2 shows a schematic diagram of a wheel speed sensor 4 
comprising a toothed wheel 5 with seven identical teeth 6. A 
sensor 7 is located at the circumference of the toothed wheel 
5 5. The sensor 7 is arranged to generate a sensor signal 
whenever a tooth 6 of the toothed wheel 5 passes the sensor 7 . 
The sensor 7 may be an optical sensor, a magnetic sensor or any 
other appropriate type of sensor which is able to detect the 
presence and the non-presence of a tooth 6 . 

10 The sensor 7 may either generate a sensor signal whenever the 
sensor 7 detects a change of its environment, i.e. whenever a 
tooth 6 of the toothed wheel 5 enters or leaves the sensor 
region, or only when a tooth 6 enters (or alternatively leaves) 
the sensor region. In the example of Fig. 2, there are in total 

15 seven sensor signals generated during one complete revolution 
of the toothed wheel 5. It is appreciated that, instead of the 
toothed wheel 5, the wheel speed sensor 4 may comprise any type 
of segmented rotary element 5 which generates a sensor signal 
for each passing sensor segment 6. Another example for such an 

20 segmented rotary element 5 is a slotted disk. The total number 
of segments is in the following denoted as L . L is not 
limited to the value chosen in the embodiment of Fig. 2 ( Z, = 7 ) 
but may be an arbitrary positive integer number. 

In more detail, the sensor 7 of the wheel speed sensor 4 
25 internally generates an internal signal with two possible 
states, high and low (e.g., high indicating a covered sensor 7 
and low indicating an uncovered sensor 7) , which in turn 
triggers the output of a clock signal delivered from a timer 
unit (not shown) , and outputs a data stream. The data stream 
30 comprises data samples in form of, for instance, a real or 
integer number t(ri) which is representative of the time 
instance of the occurrence of a corresponding internal signal. 
The time span At(n) = t{ri) — t(n -1) is defined as the duration of 
time between two successive internal signals. Thereby, n is an 
35 integer number which denotes the sample number, i.e. n = 1 
corresponds to the first sensor signal, n=2 to the second 
sensor signal, etc. 
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In Fig. 2, the solid line represents an ideal rotary element 5 
which comprises seven identical segments 6, wherein each of the 
segments 6 covers the angle a depicted in Fig. 2. The dotted 
line in Fig. 2 represents an unideal rotary element 5 in which 
5 the individual segments 6 do not have the same length but 
differ in length by an error angle 5 . These deviations from a 
nominal angle a could for example arise due to fabrication 
errors or wear during usage. In the following, the deviations 
8 from the nominal value are called imperfection errors and it 
10 is assumed that each of the segments 6 of the rotary element 
has its own characterising imperfection error S t (/ = 1,...,L). For 

instance, embodiments for estimating the imperfection errors 
are disclosed in PCT/EP02/12409 of the same applicant. The 
content of this document is incorporated into the present 
15 description by reference. In the following, a further 
embodiment for estimating the imperfection errors is described 
in more detail which is based on the embodiments disclosed in 
PCT/EP02/12409 . 

Thus, the occurrence of a sensor signal indicates that the 
20 rotary element 5 has rotated around an angle of a-2njL , in the 
ideal case of no imperfection errors, ^.nd around an angle of 
a + S t , in the realistic case with imperfection errors. From 
these sensor signals representing time instances t(ri) a 
corresponding wheel speed value coin) can be derived via the 
25 relation 

(V(n) = ^±£i (Eq. 1) 

tin)- tin -X) 

wherein a high value of coin) indicates a fast rotating wheel 
and a low value of coin) is indicative of a slowly rotating 
wheel. Besides, an estimation value for the vehicle velocity 
30 can be obtained by relating the wheel speed coin) to the 
corresponding tire radius. 
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In the following embodiments, the values t{ri) , At(n) and co(n) , 
for simplification, are all denoted as wheel speed signals and 
are considered as originating from the wheel speed sensor 4. 

For exemplification, Fig. 3 shows a diagram of wheel speeds as 
5 a function of the time, wherein the plotted wheel speeds were 
obtained during a test drive of a four-wheeled car. The diagram 
comprises four lines, each line representing one of the four 
wheels of the car. The diagram shows that during the 60 seconds 
sample period, the vehicle was driving with nearly constant 
10 velocity corresponding to a mean wheel speed of approximately 
42.3 rad/s. The diagram shows that although driving with nearly 
constant velocity the wheel speed signals are fluctuating due 
to, for example, the road roughness and the sensor 
imperfections . 

15 Fig. 4 shows, in an idealised way neglecting the influence of 
the road condition, the impact of the segment imperfections of 
a wheel speed sensor 4 on the obtained wheel speed signal a>(ri) • 
The diagram of Fig. 4 shows the wheel speed values co(n) as a 
function of the sample number n* There are 15 samples w = l,...,15 

20 shown in the diagram which correspond to three complete 
revolutions of a rotary element 5 comprising L-5 segments 6 
in total. Fig. 4 represents the case of a car 1 driving with 
exactly constant velocity v, wherein the dotted curve 
corresponds to the wheel speed signal 0){n) obtained from a 

25 wheel speed sensor 4 having an ideally segmented rotary element 
5 and the solid curve corresponds to the case of an unideal 
segmented rotary element 5 which generates a periodical 
fluctuation of the wheel speed around the average value of 5 6 
rad/s. The value of 55 rad/s of the first sample corresponds to 

30 a segment which is slightly larger than a nominal segment thus 
producing a wheel speed value which is smaller than the 
expected value of 56 rad/s. The third sample corresponds to a 
segment which exactly corresponds to a nominal segment thus 
producing the expected value of 56 rad/s. The fourth sample 

35 corresponds to a segment which is smaller than a nominal 
segment thus producing a wheel speed which is larger than the 
nominal value of 56 rad/s. The 5 th sample corresponds to the 
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last segment of the rotary element and the 6 th sample 
corresponds again to its first segment. In result, the solid 
curve of Fig. 4 shows a periodicity of five sample points which 
corresponds to a complete revolution of the rotary element 5 of 
5 the wheel speed sensor 4 . 

Below, further components of the system for estimating the road 
condition under a vehicle are explained in detail. It should 
however be noted that the subdivision of the components in 
sections and subsections has to be regarded as exemplary and 

10 not limiting. The subdivision is mainly used in order to 
increase the comprehensibility of the following embodiments. 
For the skilled person, this subdivision may also serve as a 
guideline for implementing the system. But, of course, other 
ways of structuring the system's functionality are also 

15 contemplable . Therefore, the subdivision according to the 
presented embodiments should be regarded as rather artificial 
and not as defining physical entities which can easily be 
distinguished within the final product. 

Analyser unit 

20 Fig. 5 schematically shows the components of an embodiment of 
the system for estimating the road condition. The wheel speed 
signal t(ri) obtained from the wheel speed sensor 4 is input to 
an analyser unit 8 which derives a first estimation value r(ri) 
from the received wheel speed signal t{ri) • 

25 In general, the analyser unit 8 provides an output signal (e.g. 
the first estimation value r{ri) ) which is indicative of the 
road condition under a wheel of the vehicle 1 on the basis of 
the received wheel speed signals (e.g. t{ri) or co{ri) ) of the 
associated wheel speed sensor 4 . The output signal may for 

30 example be a binary signal which indicates a rough road 
condition with a logical one (true) and a normal road condition 
with a logical zero (false) . The output signal could also be a 
real value, e.g. in the range from zero to one, whereby the 
value one indicates a maximal rough road condition, zero 

35 indicates an ideally smooth road condition and the intermediate 
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values to indicate road conditions which lie in-between these 
two extremes . 

A first embodiment of the analyser unit 8 shown in Fig. 5 
comprises a sensor imperfection estimation section 9 for 
estimating the sensor imperfections S l of the rotary element 5 
of the corresponding wheel speed sensor 4. It outputs a sensor 
imperfection signal S t which comprises sensor imperfection 

values S t , one for each segment 6 of the rotary element 5 . In a 
signal correction section 10, this sensor imperfection signal 
8 l is used to derive an imperfection-corrected sensor signal 
s(n) from the wheel speed signal t(ri) . A ground condition 
estimation section 11 then determines the first estimation 
value r(ri) of the analyser unit 8 on the basis of the 
imperfection-corrected sensor signal e{ri) . The functionality of 
the imperfection estimation section 9, the signal correction 
section 10 and the ground condition estimation section 11 is 
explained in more detail below with reference to particular 
embodiments of these sections. 

It should be noted that the above structure represents only one 
20 particular embodiment of an analyser unit 8. A second 
embodiment of the analyser unit is described with reference to 
Fig. 12 which has a different internal structure. 

Sensor imperfection estimation section 

As stated above, the sensor imperfection estimation section 9 
25 estimates the sensor imperfections 8 t of the segmented rotary 
element 5 from the wheel speed signal t{n) . 

In one embodiment of the sensor imperfection estimation section 
9, the estimated sensor imperfections S t are computed as 
weighted average values of sensor imperfection values y{n) of 
30 previous n-1 and current revolutions n of the rotary element 
5 . 
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A weighted average value may for example be obtained by a low 
pass filter which is implemented according to the following 
filter relation: 

LP : S inmoAL>x = (1 - ju)S iHnoAL>l + fiy(n) , (Eq . 2 ) 

5 with 

wherein (wmodL) + l is the number of the segment 6 of the rotary 

element 5 which corresponds to the sample number n , S nmodL is 
the estimation value of the corresponding sensor imperfection, 
10 fi is a forgetting factor of the filter, t(ri) and t{n-X) are 
consecutive values of the wheel speed signal, L is the total 
number of segments 6 of the rotary element 5 and T LAP (n) is the 
duration of a complete revolution of the rotary element 5. 

Signal correction section 

15 As stated above, the signal correction section provides an 
imperfection-corrected sensor signal e(n) based on the wheel 

speed signal t(n) and the sensor imperfection signal S l . It is 
important to note, that the imperfection-corrected sensor 
signal e(n) does not necessarily contain values which represent 

20 time instances or rotational speeds or similar quantities. It 
may also be any other artificial quantity which can 
appropriately represent an imperfection-corrected derivative of 
the wheel speed signal . 

In one embodiment, the imperfection-corrected sensor signal 
25 e{ri) is obtained from the relation 

£(n) = y(n)-S inmodL) + l (Eq. 4) 

wherein, as for the sensor imperfection estimation section 9 
(cp. above) , 
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K») = ^7TWn)-/(»-i))-^ 

wherein (wmodZ,) + l is the number of the segment 6 of the rotary 
element 5 which corresponds to the sample number n , S (nmodL)+l is 
the estimation value of the corresponding sensor imperfection, 
5 ju is a forgetting factor of the filter, t(n) and t(n-\) are 
consecutive values of the wheel speed signal, L is the total 
number of segments 6 of the rotary element 5 and T^pin) is the 
duration of a complete revolution of the rotary element 5. Of 
course, if this embodiment is implemented in combination with 
10 the embodiment of the sensor imperfection estimation section it 
is possible to use the sensor imperfection values y{n) computed 

in the sensor imperfection estimation section 9 (cf. Eq. 3) as 
input to Eq. 4. 

Ground condition estimation section and subsection 

15 As stated above, the ground condition estimation section 11 
determines the output signal of the analyser unit 8 (e.g. a 
first estimation value a,(n)) which is indicative of the road 
condition under the particular wheel of the vehicle 1 with 
which the analyser unit 8 is associated. 

20 Fig. 6 schematically shows the components of an embodiment of 
the ground condition estimation section 11. In the embodiment 
of Fig. 6, the imperfection-corrected sensor signal s(n) is 
input to a variance estimation section 12 which derives a 
variance a(n) from the imperfection-corrected sensor signal 

25 e{n) . This variance a{ri) may then be evaluated in a ground 
condition estimation subsection 13 which in turn may comprise a 
signal change determination section 14 and a decision section 
15. The signal change determination section 14 determines a 
signal change value CUSUMCounter(n) from the variance a(n) . The 

30 signal change value CUSUMCounter(ri) is input to the decision 
section 15 which outputs the first estimation value r(n) . 
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Besides, the ground condition estimation subsection 13 is not a 
necessary feature of the ground condition section 11. Fig. 7 
for example shows an embodiment of the ground condition 
estimation section 11 which solely comprises a variance 
5 estimation section 12 . 

Variance estimation section 

In general, the variance estimation section 12 computes a 
variance (here e.g. r 2 (n)) on the basis of a fluctuating input 

signal (e.g. the imperfection-corrected sensor signal £ ^ ) . 
10 There are several ways of implementing the variance estimation 
section 12 . 

In Fig. 6, the variance estimation section 12 is a subsection 
of the ground condition estimation section 11 but it may also 
be a subsection of other components (cf. the embodiment of Fig. 
15 12 in which it is a subsection of the second embodiment of the 
analyser unit 19) . 

The embodiment of the variance estimation section 12 shown in 
Fig. 8 determines a variance a(n) on the basis of the 
imperfection-corrected sensor signal e{ri) by using a low pass 

20 filter 16 (it should be noted that the term "variance" as used 
throughout the whole application does not refer to the standard 
mathematical definition but to an estimation value of the 
variance) . The low pass filter 16 may for example determine the 
variance a(n) of the imperfection-corrected sensor signal s(n) 

25 according to the following relation: 

a(n) = Var(e) = LP(e 2 )-LP(e) 2 , (Eq. 5) 

wherein LP(s) is a low pass filtered value of the imperfection- 
corrected sensor signal e{ri) , and LP(e 2 ) is a low pass filtered 
value of the square e 2 Qi) of the imperfection-corrected sensor 
30 signal s{n) . 



Here, the low pass filter 16 may be implemented according to 
the following filter relation: 
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LP : a{n + 1) = (1 - X)a{n) + Xe{n) , 



(Eq. 6) 



wherein a is an estimation value of the variance Var{s) , X is 
a forgetting factor of the filter, and s(n) is the 
imperfection-corrected sensor signal . 

5 Signal change determination section 

The signal change determination section 14 in general detects 
signal changes in an input signal (e.g. a{ri) or y(ri)) and to 
output a signal (e.g. CUSUMCounter{n)) which is indicative of 
changes in the input signal. 

10 In Fig. 6, the signal change determination section 14 is a 
subsection of a ground condition estimation subsection 13. In 
another embodiment to be described below (cf. Figs. 10 and 13), 
it is a subsection of a decision unit 18. 

In a first embodiment, the signal change determination section 
15 14 determines signal change values {CUSUMCounter{n)) according 
to the following relation: 

CUSUMCounter{n + 1) = mm(m?ot{CUSUMCounter(ri) + a(n) - Drift, 0),CounterLim it) , 

(Eq. 7) 

wherein a(n) is the variance obtained from the variance 
20 determination section, and Drift and CounterLimit are tuning 
parameters . 

Decision section 

The decision section 15 compares input values (e.g. the signal 
change values CUSUMCounter{n)) with predefined threshold values 

25 in order to derive a decision on the road condition. In 
general, the decision section 15 is optional (its input value 
already contains enough information on the road condition, its 
output signal only helps to interpret the input signal more 
easily) . For example, the decision section 15 may output a 

30 first signal indicating a rough road condition if the input 
value is higher than a threshold value, and a second signal 
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indicating a normal road condition if the input value is lower 
than the threshold value. In order to avoid fluctuations of the 
output signal when the input signal is fluctuating in the 
vicinity of the one threshold value, the results of the 
5 decision section 15 are preferably based on more than one 
threshold value. 

In the embodiment shown in Fig. 6, the decision section 15 is 
included in the ground condition estimation subsection 13. It 
may for example be designed to compare the signal change values 

io CUSUMCounter{n) from the signal change determination section 14 
with a first and a second threshold value set , reset and to 
output a current first estimation value r(n) indicative of a 
rough road condition if the signal change value CUSUMCounter(ri) 
is greater than the first threshold value set , a current first 

15 estimation value r(ri) indicative of a normal road condition if 
the signal change value CUSUMCounter(n) is lower than the second 
threshold value reset , and otherwise a current first estimation 
value r(ri) which is equal to the previous first estimation 
value r(n -1) . 

20 System for estimating the road condition under a vehicle having 
four wheels 

Fig. 9 and 10 present embodiments of a system for estimating 
the road condition under a vehicle 1 having four wheels as 
shown in Fig. 1. Each wheel of the vehicle 1 is equipped with a 
25 wheel speed sensor 4 . 

The embodiment of Fig. 9 comprises one analyser unit 8 for each 
wheel i = FL,FR,RL,RR (FL = Front -Left, FR = Front -Right, RL = 
Rear-Left, RR = Rear-Right) of the vehicle 1, wherein each 
analyser unit 8 provides a first estimation value a,(w) 
30 indicative of the ground condition under the respective wheel. 
A combination section 17 then combines the first estimation 
values a,(/?) provided from each of the analyser units 8 in 
order to obtain a combined first estimation value y{n) 
indicative of the road condition under the vehicle 1. 
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Fig. 10 shows an embodiment, in which the combination section 
17 is included in a decision unit 18 which internally post- 
processes the output value y{ri) from the combination section 17 
in order to output the first estimation value r(n) indicating 
5 the road condition under the vehicle. The decision unit 18 
further comprises a signal change determination section 14 (cf . 
description above with a{ri) replaced by y(n) ) which determines 
signal change values CUSUMCounter(n) on the basis of the 
combined output value y(n) from the combination section 17. The 
10 signal change values CUSUMCounter(n) may then be further 
processed in a decision section 15 to finally obtain the first 
estimation value r(n) . 

This embodiment can easily be adapted to any type of vehicle 
comprising an arbitrary number of sensor-equipped wheels. When 

15 a wheel speed signal t(n) is available for each wheel for 
example, then the estimation values derived thereof can be 
combined in a number of ways. Depending on the application, 
different types of tire combinations can be of interest. Some 
combinations of these are FL + RL to detect rough road left 

20 side, FR + RR to detect rough road right side or FR + FR + RL + 
RR to achieve high robustness. 

Combination section 

The combination section 17 may for example be implemented by 
computing the average value of its input signals, e.g. of the 
25 first estimation values a^ri) provided from the first analyser 
units 8 . 

Other methods of implementing the combination of the signals 
are conceivable. Alternatives are for instance networks of 
series expansion type (neural networks, radial basis function 

30 networks, fuzzy networks, etc.), min-function compared to a 
threshold, max- function compared to a threshold, average value 
compared to a threshold, or all individual signals are compared 
to a threshold and the decision is then made by voting. 
Naturally, several of the above listed alternatives can be 

35 combined. 
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System for estimating the road condition with two different 
types of analyser units 

Fig. 11 shows another embodiment of the system for estimating 
the road condition under a vehicle. It comprises two different 
5 analyser units, a first analyser unit 8 and a second analyser 
unit 19, operating on the same wheel speed signals t(ri) • 

The first analyser unit 8 is associated with the wheel speed 
sensor 4 and determines a first estimation value r x (ri) which is 
indicative of the ground condition on the basis of the wheel 

10 speed signal t{ri) received from the wheel speed sensor 4. 
Similarly to the first analyser unit 8, the second analyser 
unit 19 is associated with the wheel speed sensor 4 and 
determines a second estimation value r 2 (ri) indicative of the 
ground condition on the basis of the wheel speed signal t{ri) 

15 (respectively co(n) ) received from the wheel speed sensor 4 . 

A decision unit 18 determines a combined estimation value R(n) 
indicative of the ground condition on the basis of the first 
and second estimation values r x (n) , r 2 (n) from the first and 
second analyser units 8, 19, respectively. 

20 The first and the second analyser units 8,19 may be of a 
different type. In this case, slight differences in their 
properties can help to improve the performance of the system. 
For instance, if a first estimation value r x (n) which is output 
from the first analyser unit 8 shows weaknesses in different 

25 driving situations then a combination with a second estimation 
value r 2 {ri) which is output from the second analyser unit 19 may 

improve the detection performance. Of course, more than two 
analyser units can be combined. 

An option is to group the signals according to their source of 
30 origin, especially if the different types of signals require 
different signal processing algorithms. Due to the different 
properties of the different types of signals they are processed 
using algorithms especially adapted to this signal. Two or 
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several of the analyser units may be identical. To improve the 
algorithm even further quality measures can also be applied. 

Second analyser unit 

The second analyser unit 19 of the embodiment shown in Fig. 12 
5 comprises a band pass or high pass filter section 21 for band 
pass filtering (eg. in the range of 30-60 Hz) or high pass 
filtering the wheel speed signal a>(n) in order to remove the 
low frequency content of the wheel speed signal co(n) , such as 
vehicle acceleration. The implementation of the high pass 

10 filter may be similar to the one described in connection with 
Fig. 8. The filtering is motivated by the fact that a rough 
road, in particular a gravel road, adds (white) noise to the 
frequency spectrum of the wheel speed signal co(n) . 
Alternatively, instead of directly using the wheel speed 

15 signals co{n) already imperfection-corrected wheel speed signals 
may be used as input for the band pass or high pass filter 
section 21. The second analyser unit 19 further comprises a 
variance estimation section 12 for determining a variance value 
P(n) from the filtered wheel speed signal 3(n) , wherein the 

20 variance value f}{ri) is indicative of the ground condition under 
the respective wheel and thus can be used as a second 
estimation value r 2 (ri) which is output from the second analyser 
unit 19. The variance estimation section 19 may for example be 
similar to the one of the embodiment described in connection 

25 with Fig . 6 . 

Further embodiments of the second analyser unit 19 are 
conceivable to compute the estimation value r{ri) - For example, 
a side-wise correlation may be utilized between the front (FL 
or FR) and the rear wheel (RL and RR, respectively) on the same 
30 side of the car 1. If the vehicle moves on a rough surface, 
then the correlation at a certain velocity dependant time delay 
will be higher. An estimation value r(ri) can be obtained from 
the relations : 

R(n, k) = o)^ (w)a^ {n - k) 

r(n) = maxR(n,k) ' (Eq. 8) 
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wherein k is the sample number. A nominal value of k can be 
computed with 

Kor^al=-^TT s , (Eq. 9) 

v(w) 

where B is the distance between the front and rear axle, v{n) 
5 is the velocity of the vehicle, and T 5 is the sample period of 
v. R(n,k) can then be computed in a neighbourhood to k noYt6nal . For 

more details on correlation analysis, reference is made to 
PCT/EP03/07282 . 

Alternatively, an axle-wise correlation between the left and 

10 the right side of the car 1 may be used to determine the 

estimation value r(ri) . For a front wheel driven car the 
relation 

r(ri) = co^ (n) - (D^ {ri) - [co^ (n - 1) - oj m (n - 1)] =: (ri) - (ri) ( Eq . 10) 

may for example be used. The estimation value r(ri) is then 
15 compared to a pre-defined threshold to determine a rough road 
condition. Alternatively, the sum 

r(n)= £M",(")) (Eq. 11) 

i=FL,FR t RL,RR 

of the variance of the quantities a f (/i) defined in Eq. 10 or 
any linear combination of a subset of the four quantities can 
20 be used. In Eq. 11, Var is the variance of the quantity. 

In another alternative embodiment of the second analyser unit 
19, the analyser unit 19 monitors the highest Fourier frequency 
of the wheel speed signal according to the relation 

25 r(n) = Y J (-V k c»(k) . (Eq. 12) 

k 

The estimation value r(n) is then compared to a pre-defined 
threshold to determine a rough road condition. 



Yet another alternative embodiment of the second analyser unit 
19 can be based on the band pass filtered wheel speed signals 
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and the slip variance parameter obtained from a wheel radius 
analysis (cf . PCT/EP03/07283 ) and/or a road friction analysis. 

System for estimating the road condition under a vehicle having 
5 four wheels by means of two different types of analyser units 

Fig. 13 shows a further embodiment of the system for estimating 
the road condition under a vehicle. The embodiment is directed 
to a car 1 with four wheels j = FL,FR,RL,RR each equipped with a 
wheel speed sensor 4 . The wheel speed sensors 4 provide the 
io wheel speed signals t.(n) where i = FL,FR,RL,RR . 

One first analyser unit 8 is associated with each of the wheels 
i - FL,FR,RL,RR wherein each first analyser unit 8 provides a 

first estimation value a t (n) indicative of the ground condition 
under the respective wheel. 

15 A first combination section 17 combines the first estimation 
values provided from each of the first analyser units 8 

in order to obtain a combined first estimation value y{n) 
indicative of the road condition under the vehicle. The 
combined first estimation value y{n) is input to a signal 

20 change determination section 14 which determines signal change 
values CUSUMCounter(n) on the basis of the combined first 
estimation values y{n) according to the following relation (cf . 
above) : 

CUSUMCounter(n + 1) = mm(max(CUSUMCounter(n) + y{n) - Drift, 0),CounterLimit) , 

25 wherein Drift and CounterLimit are tuning parameters. 

One second analyser unit 19 is associated with each wheel 
i = FL,FR,RL,RR of the vehicle 1, wherein each second analyser 

unit 19 provides a second estimation value J3 f (n) indicative of 
the ground condition under the respective wheel. 

30 A second combination section 17 combines the second estimation 
values provided from each of the second analyser units 19 



21 

in order to obtain a combined second estimation value r 2 (ri) 
indicative of the road condition under the vehicle. 

An output combination section 22 finally combines the signal 
change values CUSUMCounter(ri) and the second combined estimation 

5 values r 2 (ri) in order to obtain a combined estimation value Q(«) 
indicative of the road condition under the vehicle 1. For 
instance, it may simply multiply both values CUSUMCounter{n) and 

r 2 (n) . Naturally, other signal combinations are conceivable 
(averaging, adding, etc.). The output combination section 22 

10 may be implemented similar to the first and second combination 
sections 17 and 17' as described above, in particular by a 
network of series expansion type (fuzzy or neural networks) , 
designed (trained) is such a way that it outputs a value 
between 0 and 1, with 0 representing maximum smoothness and 1 

15 representing maximum roughness. In general, all input values 
having for example values between 0 and 1 (such as a f (w) , Pi(ri) , 

y(ri) , r 2 (ri) ) may be combined with each other according to the 
above procedure . 

Optionally, a decision section 15 may be added in order to 
20 post -process the output signal Q(«) of the output combination 
section 15. An appropriate embodiment of the decision section 
15 is described above under the item "Decision section" . 

Alternative embodiment of the system for estimating the road 
condition under a vehicle having four wheels with two different 
25 types of analyser units 

Fig. 14 shows an alternative embodiment of the above system 
shown in Fig. 13. It differs from the one of Fig. 13 in that 
the signal change determination section 14 is coupled to the 
output combination section 22 instead of the first combination 
30 section 17. 

Operation results 



Fig. 15a-e show operation results of the system corresponding 
to the embodiment of Fig. 13. On the abscissas of all diagrams 
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in Fig. 15a-e are plotted the operation time interval of around 
105 minutes. On the ordinates are plotted different signal 
values obtained during the system operation. Since only the 
qualitative behaviour of the signal is relevant here, the 
magnitude of the plotted values is not defined and described in 
detail. 

In the diagram of Fig. 15a, the combined first estimation value 
y{ri) from the first combination section 17 is plotted as a 
function of the time. The diagram further shows a choice of the 
tuning parameter drift in relation to the combined first 
estimation value y(ri) . As can be seen from a shaded area in the 
diagram which represents a rough road, the first estimation 
value y{n) is larger than the tuning parameter drift on a rough 
road, and, otherwise, smaller. 

The diagram of Fig. 15b shows the signal change signal 
CUSUMCounter(ri) which is output from the signal change 
determination section 14. In principle, this signal can already 
be used to detect the road condition (for instance, if this 
signal is compared with a value CUSUMCounter{ri) = 5 as a 
threshold value) . 

The diagram of Fig. 15c shows the combined second estimation 
value r 2 (ri) which is output from the second combination section 
17. Again, this signal may already be used to determine the 
road condition . 

Fig. 15b and Fig. 15c show an interesting relation between the 
two signals CUSUMCounter{n) and r 2 (n) . They both indicate rough 
road correctly but do not incorrectly indicate rough road 
simultaneously. At 95 minutes for instance, the signal change 
signal CUSUMCounter{n) gives a strong rough road indication but 
this is not the case for the combined second estimation value 
r 2 {ri) . The opposite behaviour is present at approximately 18 
minutes . 
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The diagram of Fig. 15d shows the product Q(«) of the two 
indicators CUSUMCounter{n) and r 2 (n) as well as the two 
thresholds set and reset used in the decision section 15. In 
this diagram, a rough road is correctly indicated, whereas a 
5 rough road is not falsely indicated on a smooth road. 

The combined estimation value R(n) output from the decision 
unit 18 is shown in the diagram of Fig. 15e. Clearly, the rough 
road condition is correctly estimated in the time range from 
approximately 3 0 to 4 0 min. 

10 Computer program 

The embodiments of the computer program products with program 
code for performing the described methods include any machine- 
readable medium that is capable of storing or encoding the 
program code. The term "machine-readable medium 11 shall 

15 accordingly be taken to include, but not to be limited to, 
solid state memories, optical and magnetic storage media, and 
carrier wave signals. The program code may be machine code or 
another code which can be converted into machine code by 
compilation and/or interpretation, such as source code in a 

20 high-level programming language, such as C++, or in any other 
suitable imperative or functional programming language, or 
virtual -machine code. The computer program product may comprise 
a data carrier provided with the program code or other means 
devised to control or direct a data processing apparatus to 

25 perform the method in accordance with the description. A data 
processing apparatus running the method typically includes a 
central processing unit, data storage means and an 1/0- 
interface for signals or parameter values. 

30 Thus, a general purpose of the disclosed embodiments is to 
provide improved methods and products which enable to more 
accurately determine a rough road condition by means of wheel 
speed sensors which are in particular already existing within 
common vehicle electronic systems (antilock braking system and 

35 the like) . 
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All publications and existing systems mentioned in this 
specification are herein incorporated by reference. 

Although certain methods and products constructed in accordance 
5 with the teachings of the invention have been described herein, 
the scope of coverage of this patent is not limited thereto. On 
the contrary, this patent covers all embodiments of the 
teachings of the invention fairly falling within the scope of 
the appended claims either literally or under the doctrine of 
10 equivalents. 



